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We have performed powder inelastic neutron scattering measurements on the unconventional su¬ 
perconductor /3-FeSe (T c ~ 8K). The spectra reveal highly dispersive paramagnetic fluctuations 
emerging from the square-lattice wave vector (tt, 0) extending beyond 80 meV in energy. Measure¬ 
ments as a function of temperature at an energy of ~13meV did not show any variation from T c 
to 104 K. The results show that FeSe is close to an instability towards (tt, 0) antiferromagnetism 
characteristic of the parent phases of the high-T c iron arsenide superconductors, and that the iron 
paramagnetic moment is neither affected by the orthorhombic-to-tetragonal structural transition at 
T s ~ 90 K nor does it undergo a change in spin state over the temperature range studied. 


Iron selenide (f3- Fei +a ,Se, hereafter denoted “FeSe”) 
is structurally the simplest of the iron-based supercon¬ 
ductors but it is also one of the most intriguing. The 
superconducting transition temperature of the pure bulk 
phase is relatively low, T c « 8K pQ, but it increases to 
37 K under pressure [2] and rises above 40 K with inter¬ 
calation of alkali ions A + to form A x Fe 2 - y Se 2 [3] or by 
co-intercalation of ammonia molecules and amide ions or 
organic molecules along with A + (4 6 . Very recently, 
superconductivity was reported at temperatures as high 
as 100 K in monolayers of FeSe on SrTiOs [7, 5j- Al¬ 
though there is evidence that superconductivity at ambi¬ 
ent pressure is favored by reduction of Fe below the +2 
oxidation state and minimisation of vacancies in the FeSe 
layers 0E1, there is currently no simple explanation for 
such an extraordinary variation in T c among derivatives 
containing very similar antifluorite layers of FeSe. 

The structural and electronic ordering properties of 
FeSe differ qualitatively from those of the related iron 
pnictide compounds in two important ways. First, su¬ 
perconductivity appears in FeSe without the need for 
doping and is very sensitive to composition m- Second, 
FeSe has a tetragonal-to-orthorhombic structural transi¬ 
tion (T s ~ 90K [DE]), as in the parent phases of the iron 
pnictide superconductors, but this transition is not fol¬ 
lowed by the development of long-range magnetic order 
[1J]. The phase below T s is considered to be some form of 
electronic nematic, but opinions divide over whether the 
nematic transition is driven by orbital ordering [IMS] 
or by spin degrees of freedom pJH2Ql. 

This paper reports measurements of collective param¬ 
agnetic spin fluctuations in FeSe. Spin fluctuations are 
a prominent feature of the iron-based superconductors 
and are thought to play a significant role in the pairing 
interaction [2TH23] . In the iron arsenide superconduc¬ 
tors, spin fluctuations emerge from the same (or nearly 
so) characteristic in-plane wave vector q m = (7r,0), re¬ 
ferred to the Fe square sub-lattice, as the spin density 


wave (SDW) order of the parent phases. This magnetic 
instability is understood to be assisted by nesting of hole 
and electron Fermi surface pockets centred around the T 
and X points of the square lattice. Spin fluctuations have 
also been observed in the superconducting iron selenides, 
but the characteristic wave vector varies from system to 
system. For example, it is ( 7 r, 0) in FeTei-^Se^ (x ~ 0.5) 
[24] . (7r,7r/2) in A x Ye 2 -y ^2 (A = K, Rb, Cs) [251427] . 
and different again in Li a ,(ND 2 ) ?/ (ND 3 )i_ ?/ Fe 2 Se 2 [28] . 

Ab initio electronic structure calculations indicate that 
FeSe is close to a magnetic ordering instability with char¬ 
acteristic wave vector (tt, 0) [29H3T] . However, angle- 
resolved photoemission spectroscopy and quantum oscil¬ 
lation studies have revealed that the Fermi surface devi¬ 
ates significantly from the predictions [H1I51I32H35], and 
several models for the nematic phase predict competing 
magnetic phases with q m = ( 7 r,£), 0 < £ < 7r/2 [171420] . 
Experimental information on the magnetic ground state 
of FeSe is currently lacking, and is urgently needed to 
elucidate the nematic phase and to assess the role of spin 
fluctuations in the superconducting state. 

Here we report observations of the wave vector and en¬ 
ergy dependence of the spin fluctuations in FeSe by pow¬ 
der inelastic neutron scattering. We find collective spin 
fluctuations emerging from (it, 0) and equivalent square- 
lattice wave vectors, extending to energies greater than 
80 meV. We do not observe any significant change in the 
low energy (^10-15 meV) part of the spectrum on cross¬ 
ing the orthorhombic-to-tetragonal transition. 

A powder sample of FeSe of total mass 13.8 g was pre¬ 
pared in five separate batches of 2-3 g each. All handling 
was carried out in an argon atmosphere. Iron and sele¬ 
nium powders (5N purity) were ground together, sealed 
under vacuum in a silica glass ampoule and reacted at 
700° C for 24 h. The product of this reaction was re¬ 
ground, resealed under vacuum, annealed at 700° C for 
38 h and then cooled to 400°C and held for 6 days. The 
ampoule was then quenched in ice water and the sample 


2 



T ( K) T ( K ) d spacing (A) T(K) 


FIG. 1. (color online), (a) Magnetic susceptibility of FeSe powder. The field-cooled (FC) and zero-field cooled (ZFC) curves 
confirm the onset of superconductivity at T c ~ 8K (left). The tetragonal-to-orthorhombic structural transition at T s ~ 90K is 
signalled by a broad magnetic anomaly (right), (b) Rietveld refinement against room temperature neutron powder diffraction 
data of FeSe. Peak positions for the /3-FeSe phase are marked by vertical red ticks beneath the data. The other ticks indicate 
peak positions for Fe impurities and the vanadium sample can. (c) Temperature dependence of the orthorhombic lattice 
parameters of FeSe. The points at 150 K are the tetragonal parameters with a multiplied by y/2. The lines are visual guides. 


ground to a fine powder. The batches were found to be 
of very high phase purity by x-ray and neutron diffrac¬ 
tion, with trace amounts (< 1%) of hexagonal <a-FeSe 
and unreacted Fe as the only detectable impurities. 

Magnetisation measurements, performed with a Super¬ 
conducting Quantum Interference Device (SQUID) mag¬ 
netometer, confirmed the onset of superconductivity at 
T c ~ 8 K in each of the five batches. An example of field- 
cooled and zero-field-cooled data is shown in Fig. [lja). 
Measurements retaken after the neutron scattering ex¬ 
periment confirmed that the sample did not deteriorate. 
The right-hand panel of Fig.JIJa) shows a broad magnetic 
anomaly at the structural transition T s 90 K consistent 
with previous data on FeSe powders pQ. 

For a detailed structural analysis, we performed high 
resolution neutron powder diffraction on the HRPD in¬ 
strument at the ISIS Facility. Measurements were made 
at temperatures between 10 K and room temperature. 
Figure [ljb) shows data collected at room temperature 
together with a Rietveld fit. The temperature depen¬ 
dence of the lattice parameters obtained from the refine¬ 
ments are shown in Fig. [ljc). The continuous tetrago¬ 
nal ( P4/nmm ) to orthorhombic ( Cmma ) transition at 
T s ~ 90 K is consistent with earlier results mmm- 
The orthorhombic distortion (b — a)/a approaches 0.5% 
at 10 K. Refinement of the composition Fei +a; Se against 
data above and below T s yielded x = 0.01(1), i.e. with 
interstitial Fe sites between the stoichiometric FeSe lay¬ 
ers occupied at the 1% level, a finding consistent with a 
previous report correlating composition with t c m- 

Inelastic neutron scattering was performed on the 
chopper spectrometer MERLIN at the ISIS Facility [37]. 
The powder sample was loaded into aluminium foil pack¬ 
ets and placed in an aluminium can in annular geome¬ 
try. The can was attached to a closed-cycle refrigera¬ 
tor. Neutron spectra were recorded with incident ener¬ 


gies of E[ = 34, 50 and 100 meV at temperatures from 8 
to 104 K. The spectra were normalised to the incoherent 
scattering from a standard vanadium sample measured 
with the same incident energies, enabling us to present 
the data in absolute units of mb sr -1 meV _1 f .u. -1 (where 
f.u. refers to one formula unit of FeSe). 

Figure [2|a) shows an intensity map of part of the 
Ei = 100 meV spectrum measured at 8 K on MERLIN. 
The spectrum is dominated by scattering from phonons 
for energies E below the phonon cut-off at 40meV [38] . 
Above 40 meV, there is a broad vertical column of scat¬ 
tering centred on the wave vector Q = 2.6 A -1 , and a 
weaker column centred on 3.5 A -1 . Figure (2^b) is a sim¬ 
ilar intensity map measured with E\ = 34 meV to probe 
the low (Q, E) part of the spectrum. Phonon scattering 
dominates in this regime, but there is a window between 
10 and 15 meV in which the phonon signal is small, and 
a vertical column of weak scattering can be seen centred 
near Q = 1.2 A -1 . Such scattering columns are observed 
in neutron powder spectra of other iron-based supercon¬ 
ductors and have been confirmed to arise from strongly 
dispersive spin fluctuations j39l - f42j . 

The magnetic signals identified in the intensity maps 
can be seen in more detail in the Q cuts made at fixed 
average energy shown in Fig. Ib). The cuts contain 
peaks centered on Q = 1.2, 2.6 and 3.5 A -1 , and there 
are additional weak signals near Q = 4.5 A -1 . The series 
of magnetic peaks can be indexed as orders of the square 
lattice wave vector (7r,0), see Fig. |3|a). In reality, the 
magnetic signal will extend in the out-of-plane direction, 
either as a diffuse rod of scattering if the correlations 
are quasi-two-dimensional or as a series of peaks if there 
are strong inter-layer correlations. Simulations of such 
types of out-of-plane scattering show that after powder 
averaging the peaks have a tail on the high Q side but 
the maxima shift by only a small amount (< 0.06 A -1 ) 
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FIG. 2. (color online). Neutron powder spectra of FeSe 
obtained on MERLIN at T = 8 K ~ T c . The vertical dashed 
lines show the 2D wave vector ( 7 r, 0) and equivalent positions, 
(a) High-energy part of the spectrum recorded with an inci¬ 
dent energy E\ — 100 meV. The vertical bands of scattering 
above the phonon cut-off at 40meV are caused by steeply 
dispersing cooperative paramagnetic fluctuations, (b) Low 
energy part of the spectrum from the data measured with 
Ei = 34meV. Magnetic scattering is visible in the energy 
window 10-15 meV where phonon scattering is weak. 


from the ideal two-dimensional wave vectors. 

Although FeSe does not order magnetically, our results 
show that it has a strong magnetic response at (7r, 0) and 
equivalent positions which characterise the in-plane SDW 
order found in the parent phases of the iron arsenide su¬ 
perconductors. To quantify the analysis, we compare the 
data to a phenomenological model for the low energy re¬ 
sponse of a two-dimensional (2D) antiferromagnetically- 
correlated paramagnet. The model has been used previ¬ 
ously to describe the low energy part of the spectrum of 
superconducting Ba(Fei_ a ,Co a ,) 2 As 2 [43]. The neutron 
scattering cross section may be written 


d 2 cr 

dfMF/f 


^S(Q,E), 


( 1 ) 


where 5(Q,£'), the magnetic response function, is the 
quantity presented here. For an isotropic paramagnet, 


= (£) 2 l-exp(-^) ^ 3<Q)X, ' (q - E) - 

( 2 ) 

where (yro/2) 2 = 72.7mb, /3 = l/ZcpT, /(Q) is the mag¬ 
netic form factor, and x"(q, E) is the absorptive part of 
the generalized susceptibility. The low-energy magnetic 
excitations are envisaged as damped spin waves with a 



FIG. 3. (color online), (a) Reciprocal lattice of the Fe 
square lattice. The (red) filled circles mark the (77 0)-type 
wave vectors and the (green) filled diamonds the (77 7 t /4)- 
type positions, including the equivalent 90° domains. The 
full and dotted circles show the effect of powder averaging, 
and the table lists the corresponding values of Q = |q m |. (b) 
Constant-energy cuts through the data in Fig. [ 2 ] for three dif¬ 
ferent energy bands as indicated. The upper two cuts are off¬ 
set vertically. The (red) symbols are the data and the (blue) 
lines are fits with a damped harmonic oscillator model for 
spin-waves dispersing anisotropically from q m = ( 7 r, 0 ). (c) 
The same as in (b) but with q m = ( 7 r, 7 r/ 4 ) and an isotropic 
dispersion. 


linear dispersion, and we use a harmonic oscillator model 

2E 2 TE 

X (q, E) oc ^2 _ £ 2)2 _|_ 4 p 2£2 > ( 3 ) 

in which E q = h[(v\\q \\) 2 + (v±q±) 2 ] 1 / 2 is an anisotropic 
dispersion with velocities v\\ and v± in the longitudinal 
and transverse directions relative to q m = (7r,0), T = 
7 E is the inverse lifetime, and q is the spin-wave wave 
vector. x"(q, E) does not vary with q z , and is repeated 
in 2D momentum space with the periodicity of the 2D 
magnetic wave vector q m . We fitted the model to the 
constant-energy cuts allowing 7, v \\, v±, an intensity scale 
factor and a flat background to vary. The experimental 
Q resolution was included. 

The parameters determined from the fit are 7 = 0.13T 
0.06, t7|| = 460zb 120 meVA and v± = 150±20meVA. The 
anisotropic velocity obtained from this analysis is statis¬ 
tically significant. Spectra simulated with the best-fit 
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FIG. 4. (color online). Temperature dependence of the mag¬ 
netic scattering at (77 0) ~ 1.2 A -1 averaged over the energy 
range 11 — 14meV. The solid lines are fits to a Gaussian func¬ 
tion on a linear background (dotted). The upper three scans 
are offset vertically by 0.25, 0.5 and 0.75 units, respectively. 

parameters are shown in Fig. [3jb) [44]. The simulations 
match the peak at Q = 1.2 A _T and closely reproduce the 
observed dispersion of the signals centered near 2.6, 3.5 
and 4.5 A -1 . The model parameters are similar to those 
found for Ba(Fei_ a ,Co a ,) 2 As 2 : 7 = 0.15, v\\ = 580meVA, 
v± = 230meVA [43]. We also considered a model for 
purely diffusive spin dynamics [45]. The diffusive model 
also fits the data satisfactorily and leads to the same con¬ 
clusions as the damped spin wave model [44] . 

Despite the limitations inherent in powder-averaging, 
the success of the model in accounting for features in the 
data over several Brillouin zones places a tight constraint 
on the wave vector q m that describes the dominant mode 
of paramagnon excitations in FeSe. As a test of this, 
we carried out fits with the damped spin-wave model 
modified to have q m = ( 7 r, £) and an isotropic dispersion. 
The fits with this model were in all cases inferior to those 
with q m = ( 7 r, 0) and an anisotropic dispersion. The best 
agreement was achieved with £ « 7t/4 and is shown in 
Fig. § c), but the leading edge of the fitted signal near 
2.6 A -1 is at too low Q compared with the ( 77 , 0) model 
reflecting the difference between the magnitude of the 
wave vector (7t,27t), Q = 2.64 A -1 , and ( 77 , 77t/4), Q = 
2.38 A -1 — see Fig. |3[a) and [44] . 

In this experiment we were unable to cool the sam¬ 
ple below 8K, and so did not study the magnetic sig¬ 
nal in the superconducting state at low energies where 
a spin resonance could be expected. Instead, we in¬ 
vestigated the influence of the structural transition on 
the magnetic response by performing runs with E[ = 
50meV at temperatures of 8, 17, 67 and 104 K. Fig¬ 
ure [4] shows Q cuts through the ( 77 , 0) position at each 
temperature. The data are averaged over the energy 


interval from 11 to 14meV to stay within the win¬ 
dow where phonon scattering is weak. The magnetic 
peaks show very little variation with temperature. To 
quantify this, we fitted a Gaussian function on a lin¬ 
ear background to each cut. To within the fitting er¬ 
ror the integrated intensity remains constant at 0.10 db 
0.01 mbsr -1 meV -1 A -1 f.u. -1 , which compares with the 
value 0.08 d= 0.01 mbsr -1 meV -1 A -1 f.u. -1 found at the 
same energy for LiFeAs at T = 20 K > T c [42]. This 
shows that the spin fluctuations in FeSe have a similar 
strength to those in other Fe-based superconductors. 

The fact that the magnetic response shows very little 
or no change on crossing the structural phase transition 
implies that the structural transition is not driven by 
magnetic fluctuations at the frequencies probed in our 
experiment. Further, the lack of any change over the en¬ 
tire temperature range studied implies that the param¬ 
agnetic moment is constant below 104 K, in contrast with 
the notion of a gradual spin-state transition proposed to 
explain thermally-induced phonon anomalies observed in 
Raman spectra [46] . 

This study establishes that the collective spin fluctua¬ 
tions in FeSe share many similarities with those in the 
high-T c Fe arsenide superconductors, including a very 
steep dispersion and a low frequency response that is 
strongest at or very close to the square lattice wave vector 
(77 0). We find no direct evidence for competing mag¬ 
netic orders, although the highly anisotropic spin-wave 
velocity implies a greater tendency for transverse spin 
fluctuations. If spin fluctuations are important for the 
pairing mechanism in Fe-based superconductors then our 
results show that the ingredients for high-T c are present 
in FeSe, and something other than conventional magnetic 
dipole fluctuations must compete with superconductivity. 
Several different nematic degrees of freedom that could 
suppress superconductivity have been discussed recently 
H3H20], and experiments to search for possible orbital 
and spin nematic order parameters compatible with ( 7 r, 0) 
spin fluctuations will be an important next step. 

Note added : Recently, an eprint appeared reporting 
neutron scattering measurements of the low-energy re¬ 
sponse in FeSe m- 
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Damped harmonic oscillator model 

The damped harmonic oscillator (DHO) model is de¬ 
fined in Eqs. 0 0 of the main article. It is a phe¬ 
nomenological model for damped spin waves. The model 
has four variable parameters: the longitudinal and trans¬ 
verse spin-wave velocities v\\ and v± , the damping param¬ 
eter 7 , and an intensity scale factor. This model fits the 
experimental spectrum well, as demonstrated via the Q 
cuts shown in Fig.J^b) of the main article. The fit shows 
that the spin-wave velocity is anisotropic with v\\ ~ 3v±. 
Figure SI (left panels) shows intensity maps for the best- 
fit DHO model simulated over the same range of energy 
and wave vector covered by the measurements in Fig. [2] 
of the main article. 

We also consider the possibility that the spin fluctu¬ 
ations are associated with a characteristic wave vector 
q m = (7r, £), as indicated by some theoretical models. To 
perform a simple test we modified the DHO model to 
have q m = ( 7 r,£) and fitted the model to the data for 
0 < £ < 7t/2, constraining for simplicity the spin-wave 
velocity to be isotropic. Figure [3jc) in the main article 
shows the best fit that could be achieved with this mod¬ 
ified DHO model, which was with £ ~ 7t/4. The fit is 
not as good as the DHO model with anisotropic velocity 
and q m = (77 0). The discrepancies are most noticeable 
for the signal near 2.6 A -1 at low energies. The reason 
for this is illustrated in Fig. J3^a) of the main article. 
The unmodified DHO model converges at low energies 
to q m = ( 7 r, 0) and equivalent points in 2D momentum 
space, whereas the modified DHO model dispersion con¬ 
verges to two points at q m = ( 77 , ±7t/4), and equivalent 
points. Therefore, in a powder-averaged spectrum at low 
energies the anisotropic model will show a rise in inten¬ 
sity at Q = 2.64 A -1 corresponding to q m = (7T, 27r), 
whereas the isotropic model will show a rise in intensity 
at Q = 2.38 A -1 corresponding to q m = ( 77 , 77t/4). To 
illustrate this effect more directly we show in Fig. S2 the 
distribution in intensity in 2D momentum space for the 
two models calculated for the same energies as the Q cuts 
in Fig. [3j 



Q ( A ' 1 ) Q ( A ' 1 ) 


Fig. SI. Intensity maps of the best-fit DHO model (left pan¬ 
els) and diffusive model (right panels). The spectra cover the 
same range of energy and wave vectors as in Fig. [2] of the 
main article. 

Diffusive model for correlated paramagnet 

Strictly speaking, the DHO spin-wave model applies to 
systems with short-range magnetic order. The justifica¬ 
tion for applying it to FeSe is that there are strong mag¬ 
netic correlations, which means that the spins will appear 
ordered over sufficiently short distances and timescales. 
Therefore, above some crossover energy scale the spec¬ 
trum is expected to be very similar to that of a system 
whose spins are ordered with propagation vector q m and 
whose excitations are propagating damped spin waves. 

An alternative approach is to employ a description that 
does not imply any magnetic order in the ground state. 
Such a phenomenological model was used by Tucker et 
al. to analyse neutron scattering data on superconduct¬ 
ing Ba(Fei_ x Co ;E ) 2 As 2 with x = 0.047 [45]. This model 
is designed to capture the diffusive nature of the spin dy¬ 
namics inherent in a paramagnet while at the same time 
allowing spatial and temporal magnetic correlations. The 
scattering intensity is proportional to the imaginary part 
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Fig. S2. Intensity maps of the DHO model in 2D momentum 
space at different energies. The magnetic wave vector is q m = 
(tt, 0) in the left panels (best-fit model), and q m = (7r,±7r/4) 
in the right panels. 


In Eqs. (E])-([6]), q = (q\\,q±) is the reduced wave vector 
measurea relative to a magnetic wave vector, i.e. q = 0 
is at q m . The components q\\ and q± are parallel and 
perpendicular to the vector joining q m to the nearest 
reciprocal lattice vector, e.g. for q m = (7r,0), q\\ and q± 
are components parallel and perpendicular to q m . At 
low energies this model describes an elliptical intensity 
distribution centred on q = 0, with q\\ and q± as the 
principal axes of the ellipse. x"(q, E) is repeated in 2D 
momentum space at each of the 2D wave vectors q m . 

The best fit of the diffusive model to the data was ob¬ 
tained with parameters = 7 =b 2 A, ^ = 1.5 =b 0.8 A, 
T|| = 0 meV, r ± = 7 ± 5 meV and a = 7 db 5 A. The pa¬ 
rameter T|| was not well controlled and had little effect 
on the fit. The large difference between and £_l means 
that the intensity distribution is highly anisotropic, con¬ 
sistent with what was found with the DHO model. 

Spectra simulated with the best-fit parameters of the 
diffusive model are shown in Fig. S3. The simulated Q 
cuts match the data well when q m = (7r, 0) (left panel), 
but less well with q m = (7r,±7r/4) and an isotropic in¬ 
tensity distribution (right panel). 

The main conclusions are, firstly, that the phenomeno¬ 
logical diffusive model provides a good description of the 
magnetic dynamics of FeSe, and second, that the analy¬ 
sis with the diffusive model reinforces the findings from 
the DHO spin-wave model analysis that FeSe is close to 
an instability towards (n, 0) antiferromagnetism and has 
highly anisotropic magnetic correlations. 


of the generalized susceptibility, see Eqs. ([l]) ([ 2 ]), which 
in the diffusive model is given by 


x"(q ,E) <x 


r q E 

Et + Tld+Qq 2 ) 2 " 


(4) 


where the momentum-dependent correlation lengths and 
relaxation rates are 


&z 2 «£f«f + &zl (5) 

r q = r + ^(rf< 7 f + ri (/ i). (6) 

The model has six variable parameters including an in¬ 
tensity scale factor. The parameters used here are consis¬ 
tent with those defined by Tucker et al. [45{. Specifically, 
the parameters £, T, 77 ^ and 77 p in Tucker et al. are given 
by 



— C 2 (i + 


(7) 

Fig. S3. Constant-energy cuts through the data together with 
fits with the diffusive model, (a) q m = (7r, 0), (b) q m = 

rf =r 2 (i + ?7r), 

r i = r 2 (l-7?r)- 

(8) 

(7t,7t/ 4) and an isotropic dispersion. 



















